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Abstract. - We show how a recent proposal to obtain the distribution of conductances 
in three dimensions (3D) from a generalized Fokker-Planck equation for the joint probability 
distribution of the transmission eigenvalues can be implemented for all strengths of disorder 
by numerically evaluating certain correlations of transfer matrices. We then use this method 
to obtain analytically, for the first time, the 3D conductance distribution in the insulating 
regime and provide a simple understanding of why it differs qualitatively from the log-normal 
distribution of a quasi one-dimensional wire. 



It has been shown recently that for a quasi one-dimensional (QID) disordered quantum 
wire, the distribution of conductances P[g) (at zero temperature, in the absence of electron- 
electron interaction) has many surprising features arising from large mesoscopic fluctuations. 
These include a highly asymmetric 'one-sided' log-normal distribution at intermediate disorder 
between the metallic and insulating limits [1,2], and a singularity in the distribution near the 
dimensionless conductance 17 ~ 1 in the insulating regime [3]. Although there is no phase 
transition in QID conductors, numerical studies support the conjecture that some of these 
features persist in higher dimensions as well [4, 5] , which may have important consequences 
for the Anderson metal-insulator transition in three dimensions (3D) [6]. However, while a 
systematic method has recently been developed to obtain analytically the full P{g) in QID [7], 
no such method is yet available in higher dimensions. Indeed, it has not been possible so far 
to study analytically even the simpler case of the distribution P{g) for a 3D insulator [8,9]. 

A phenomenological generalization of the QID method to study the distribution of con- 
ductances in higher dimensions has been proposed recently [10]. For a given dimensionality, 
the generalization involves an unknown matrix K to be determined from numerical studies of 
the properties of transmission matrices for conductors of various strengths of disorder. Within 
a set of well-defined approximations, this method is applicable in principle for all strengths of 
disorder in all dimensions. In the present work we first of all test these approximations and 
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show that they are vahd in 3D. We then study the matrix K in 3D and show that it can be 
modeled reasonably well by a single parameter in both the metallic and insulating limits, and 
that it contains information about the critical regime. We find that while the value of the 
parameter in the metallic regime in 3D remains the same as in QID, it changes dramatically 
across the transition and in the insulating regime. This implies that while the 3D metal is 
similar to a QID metal, the localization in 3D is qualitatively different from the localization 
in QID, discussed in Refs. [1-3]. We use the value of the parameter in the insulating limit 
to obtain analytically, for the first time, the full P{g) of a 3D insulator. We show e.g. that 
the density of Lyapunov exponents in 3D is qualitatively different from that in QID; while 
the conductance in the insulating regime is dominated by the smallest Lyapunov exponent in 
QID, this is no longer true for 3D. The resulting 3D P{g) is not log-normal even in the deeply 
insulating regime. Our results are in agreement with direct numerical studies. 

The distribution of conductances in QID was obtained within the transfer matrix approach. 
In this approach, a conductor of length Lz and width L x i is placed between two perfect leads; 
the scattering states at the Fermi energy then define N (x L'^ channels. The 2N x 2N transfer 
matrix M relates the flux amplitudes on the right of the system to those on the left [11]. Flux 
conservation and time reversal symmetry (we consider the case of unbroken time reversal 
symmetry only) restricts the number of independent parameters of M to N{2N + 1) and M 
can be written in general as [12] 

" - {I {"W Jh) {I :■) ■ « 

where u, v are N x N unitary matrices, and A is a diagonal matrix, with positive elements 
Xi,i = 1,2, ...N. The restriction of the approach to QID arises from the "isotropy" approx- 
imation, that the distribution pi^^{M) is independent of the matrices u and v. In [10], this 
restriction was lifted in favor of an unknown phenomenological matrix 

TV 

Kab^ikab)^; kab = ^\Vaa\'^\vba\'^, (2) 
a 

where the angular bracket represents an ensemble average [13]. In terms of this matrix, the 
distribution (M) satisfies an evolution equation given by [10] 

J^l[\X,-X,\-y^^; lab^-p^. (3) 

a<b 

In QID, the matrix K reduces to Kab = ''n+i ; which implies ^ab = 1, and one recovers 
the well known DMPK equation with the symmetry parameter /3 = 1 [12]. In 3D, K is not 
known analytically. We obtain it from direct numerical studies of a tight binding Anderson 
model defined by the Hamiltonian 

n = W^Encicn + ^ tnn'c\cn'. (4) 
n [nil'] 

In n = {xyz) counts sites on the lattice of the size Lx Lx L^, and e„ are random energies, 
uniformly distributed in the interval [—5, 5]- The parameter W measures the strength of the 
disorder. The Fermi energy is chosen as Ep = 0.01. The hopping term t„„/ between the 
nearest-neighbor sites nn' is unity for hopping along the z direction and tnn' — t for hopping 



dpL, (A) _ j_ 

diL,/l) ~ J^dXa 



Aa(l + Xa)KaaJ 



dp 

dXa 
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Fig. 1 - Left: the distributions of normalized parameters fcn (open symbols) and fci2 (full symbols) 
in the metallic regime {W = 4). Two right figures show the same for the insulator [W — 20). Inset 
of right figure shows Ku and K12 of systems 6 x 6 x as a function of L^. 

in the x and y directions. To avoid closed channels existing in perfect leads, we use t = 0.4. 
Then the model exhibits a metal- insulator transition at Wc « 9. In order to construct 
a simple model for the matrix K, we use the method of [14] to calculate numerically the 
matrix TT'^ (T is the transmission matrix). Using Eq. (Q, we have TT^ — v*{l + X)~^v. 
Diagonalization of TT^ gives us both A and all elements of the matrix v. 

There are two major assumptions made in [10] in deriving Eq. (i) the elements kab 
can be replaced by their mean values Kab and (ii) the L^-dependence of Kab is negligible. 
To test these assumptions, we analyzed the probability distributions P {kab/ Kab) in the in- 
sulating regime and compared them with those for metals (figure 0. In the metallic regime 
all distributions seem to be self-averaging (they become narrower when L increases), with 
sharp maximum at 1. In the insulating regime P{ki2/ K12) has a very sharp distribution 
too, but possesses a long exponential tail to values fci2 ^ K12. P{kii/Kii) is broader, with 
variance of the same order as the mean. This means that while assumption (i) remains valid 
to leading order, fluctuations in the diagonal elements kaa in the insulating regime can be 
important if the final results are sensitive to the exact values of the parameters. We will 
therefore concentrate on the qualitative features of the distribution of conductances, which 
are insensitive to those fluctuations. Figure Q also shows that Kab depend slightly on the 
ratio Lz/L and reach L^-independent limiting values when LyjL — > 00. This is qualitatively 
consistent with assumption (ii) and agrees with previous numerical analysis of parameters xi 
(Ai = sinh'^Xi) [15,16]. We therefore conclude that at least to leading approximation, the 
generalized DMPK equation Q can be used in 3D at all disorder strengths. 

We now construct a model for the matrix K in the insulating regime. Figure |21 shows Kaa 
as a function of a and ^ab for different (a, b) in this regime. While there is some dependence 
on the indices for both K and 7, the dependence is weak compared to the dependence of 
the matrix elements on disorder in the L — > c» limit. As only few channels contribute to the 
conductance in the insulating regime, we will ignore this index dependence and use Kaa ~ Ku 
and 7a6 ~ 712 [17]. FigureOlis consistent with Ku oc 1/L™, where m = 2, 1,0 in the metallic, 
critical and insulating limits, respectively, in agreement with [18]. Figure 13 also shows 712 
as a function of disorder. It is unity in the weakly disordered metallic regime, showing the 
same behavior as in the QID limit. This is true for all 7a6, which implies that a 3D metal 
is very similar to a QID metal. However, as the disorder is increased for given L, the value 
of 712 starts to decrease. Several plots with different L values show that in the insulating 
limit 712 cx 1/L (right figure |2Jl. Thus our data show that the behavior of strongly disordered 
insulators differs from that of QID systems [8]. 

While modeling the full 7ab at the critical point needs more careful numerical studies, the 
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Fig. 2 - Left: a-dependence of Kaa for W — 30. Right: 71;, for two values of L and (in inset) ^aa+i 
for i = 1, 2 and 3 (W = 30, L = 10). 



insulating limit is simpler and provides a test case for the generalized DMPK equation (3). 
It predicts that the logarithmic interaction between the transmission eigenvalues Xa vanishes 
as 1/L in the insulating limit. In this paper we will test this prediction by evaluating the 
full distribution of conductances in the insulating limit for a 3D conductor as described by 
Eq. ||2J), using the simple approximate model for K suggested by our numerical studies [19], 
namely 



and 7ab « 712 ~ e/2i- 



(5) 



We can now follow [20] and map the problem onto a Schrodinger equation in imaginary time. 
The corresponding Hamiltonian contains an interaction term with strength proportional to 
712(712 — 2) which vanishes not only for the QID unitary case 712 =(3 = 2 but also in the limit 
712 — > 0. In the present case, 712 is of order ^/2L <C 1, and we can neglect the interaction 
and use the (3 = 2 solution of [20] . For Aq = sinh^ Xa the distribution is then given by [1] 



N 

Pig) = y n dxae-"5{g - ^ h{xa))\ h{x) = sech^ 



(6) 



Here the (5-function represents the Landauer formula for the conductance of a multichannel 
system [21], and the 'Hamiltonian' H is given, in the insulating limit Xa S> 1, by [20] 



N 



a>b 



— In I sinh^ Xa — sinh^ xij\ 



712 



In! 



N 

E 

a=l 



— In sinh 2a; „ 
2 



In; 



Txt 



(7) 

where F = l/LzKn « 2^1 ^ 1. The ratio r/712 ~ '^L/L^ depends on the geometry, 
but is independent of disorder, leaving a single independent parameter F that determines 
the strength of disorder. We will use F/712 = 4 appropriate for a cubic system. Note 
that fluctuations in kaa, ignored in the model, would make the numerical factors in F and 
712 inaccurate, but will not change either the length or the disorder dependence of these 
parameters. 

The replacement of /? = 2 in [20] by 712 ^ in Eq. ^ has the consequence that while 
all {Xa) ^ 1 in the insulating regime, the difference s = {xa+i — Xa) is not of the same order 
as {xa)- For example, if we keep only the first two levels, the saddle-point solutions for xi 
and X2 give (xi) ^ Lz/S, and {x2 — Xi) ^ {xi). We therefore do not assume that X2 '3> Xi. 
However, we do make the simplifying approximation that In | sinh^ Xa — sinh^ X})\ ~ Insinh^ Xa 
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Fig. 3 - Left: LKn for various L as a function of disorder. Note the common crossing point for 
~ 9, showing Kn oc 1/L at the critical point. Inset shows the L-dependence of LKn for three 
values of disorder W. Solid lines are power fits LKu oc with m = —1.3,0.06 and 0.8 for for 
W — 2 (metal), W — 9 (critical point) and W — 30 (insulator), respectively. The linear fit for 
W — 30 uses Kii — 0.37 + 1.2/L. Middle: 712 as a function of disorder for various system sizes. In 
the metallic regime, 712 ~ 1. At the critical point, 712 ~ 0.25. In the localized regime, 712 decreases 
with increasing W . Inset shows 712 (i) for W = 2,9 and 30. Solid line is a fit 712 — 0.015 + 0.56/L. 
Right figure confirms that 712 oc 1/L in the localized regime. Limiting values of 712 (L — > 00) are 
given in the legend. 

and In \x1 — « Inx^ for a > 2. Eq. iQ then becomes 

N 

H^Hi+Y,[V{xa)-Ji2{a-2)f{xa)], (8) 

a=2 

where 

Hi = — In — x^l + Fa;^ lnsinh2xi — Inxi, 

1 2 

V{x) = Fa;^ — — lnsinh2x — Ina; — 712 Insinhx; f{x)—lnsmhx-\ In a: (9) 

2 712 

Following ref [1], we separate out the lowest level xi and treat the rest as a continuum 
with density a{x) beginning at a point X2 > xi. The corresponding saddle point free energy 
Fsp{xi,X2't g) has the form 

Fs,{xui2) ^H,-— j^[V"{x) - f^lk'^x}], (10) 

^712 Jx2 J \^) 

where /xi = —V'{x2)/h'{x2) and primes denote x-derivatives. Eq @ can then be rewritten as 
P{\ng)(xg dS:2e~''"'^'=^'^^'S') e'^^^^-'''^ (11) 

where the integration over is eliminated by a constraint arising from the minimization of 
the free energy: 

XI =cosh-i[l/V^^; 30 = / ^y^M{x)[V'{x)-\-^llh\x)]. (12) 

712 Jx2 J 

The lower limit X2m\n is the larger of the additional constraints imposed by the conditions 
c(i2) > and X2> Xi> 0, Xi real. 
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Fig. 4 - Left: Conductance distribution for 3D insulators obtained from direct numerical simulation 
for W = 20, L — 18 (circles), and from Eq. IIIH for F — 0.054 (solid line). Both have the same 
mean value {Ing) = —12.6. Dashed and dotted lines show Eq. 1131 with i2min = l/2r + 5/8 and 
£2min = l/F, respectively, with F = 0.054. Right: Comparison of var(lng) and skewness as a function 
of (Ing) obtained from numerical simulations for various W and L (full symbols) and from present 
model Eq. lllll (open symbols). 

Let us consider first a simple approximate solution of Eq. (|11|) , which is dominated by the 
lower limit of the integral. To a good approximation, go is negligible compared to g in the 
insulating limit, and xi « i ln{4:/g). The condition a{x2) > gives a;2min ~ (1 + F + 7i2)/2F 
and hence Fgp « Hi. This immediately leads to 

P(lng)cx(4iLin-w')e-5^^+«)^, u = Hg/4) (13) 

valid for \u\ < 25:2min- Figure 0] shows Eq. H13|l compared with the results from direct inte- 
gration of Eq. 1)11(1 , both compared with numerical results based on Eq. . For the analytic 
curves, we chose F = 0.054 to have the same < Ing > from Eq. 1(11(1 as in the numerical case. 
Note that using the QID result 712 — 1 gives X2min ~ l/T, leading to a log-normal distri- 
bution (sec dotted line in Figure 0J. A saddle point analysis of Eq. 1(13(1 allows us to obtain 
(Ing) w l/F — ■\/2/F and var(ln5')« l/F, to be compared with the corresponding QID results 
(Ing) « l/F and var(lng)« 2/F. Moreover the skewness ((Ing— (lng))^)/[((lng — (Ing))^]^/^ 
saturates to a finite value of order unity independent of F for F ^ 0, showing the absence of 
log-normal distribution for 3D insulators even at very strong disorder. Right figure 0] shows 
variance and skewness calculated from direct integration of Eq. I(ll() compared to numerical 
results, consistent with saddle point results from Eq. ((13(1 . Quantitative differences between 
Eq. 1(11(1 and numerical results are due to our simplified model Eq. ((SJ, which still overestimates 
the strength of the interaction for higher channels. 

It is instructive to analyze the eigenvalue spectrum in terms of the density (t{x). In the 
insulating phase we find a{x) ~ 2F/712 — SL/L^ for x :s> 2/712. This corresponds to a 
uniform average spacing s — {xa+i — Xg) of eigenvalues of order unity (L ~ L^), compared 
to the uniform spacing s ~ L^/^ in QID. In contrast, 3D metals are similar to QID metals 
having uniform (j{x) extending down to a; = and s ~ Lzjl? . The opening of a gap in the 
spectrum of Lyapunov exponents Vn = {xn) / Lz ~ 1/C may be considered as the signature of 
the Anderson transition. 

We conclude that Eq. ^ indeed describes a 3D insulator. The solution of Eq. |j2Jl, given 
in Eq. ((11(1 proves the non-trivial asymmetry in the distribution P(lng), which is qualitatively 
distinct from a QID insulator. Although such a distribution has been known numerically, our 
method gives for the first time a simple theoretical understanding of the entire distribution. 
This opens up the possibility to study analytically in more detail the insulating regime as 
well as the Anderson transition in 3D in terms of the distribution of conductances, providing 
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an opportunity to investigate at least qualitatively the nature of a quantum phase transition 
in the presence of large mesoscopic fluctuations. Numerical work is underway to construct a 
model for the matrix K at the critical point. 
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